The glass transition and the replica symmetry breaking in vortex matter 
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s ! . We quantitatively describe the competition between the thermal fluctuations and the disorder 

using the Ginzburg - Landau approach. Flux line lattice in type II superconductors undergoes 
a transition into three "disordered" phases: vortex liquid (not pinned), homogeneous vortex glass 
(pinned) and crystalline Bragg glass (pinned) due to both thermal fluctuations and random quenched 
disorder. We show that disordered Ginzburg - Landau model (valid not very far from H C 2) in which 
only the coefficient of a term quadratic in order parameter tp is random first considered by Dorsey, 
Fisher and Huang leads to a state with nonzero Edwards - Anderson order parameter, but this 
state is still replica symmetric. However when the coefficient of the quartic term \i/j\ 4 in GL free 
energy also has a random component, replica symmetry breaking effects appear. The location of 
the glass transition line in 3D materials is determined and compared to experiments. The line is 
clearly different from both the melting line and the second peak line describing the translational 
5— I ' and rotational symmetry breaking at high and low temperatures respectively. The phase diagram 

Qh| is therefore separated by two lines into four phases mentioned above. 

^ ■ PACS numbers: PACS numbers: 74.20.De,74.60.-w,74.25.Ha,74.25.Dw 

£h ■ I. INTRODUCTION 

-a ' 

In any superconductor there are impurities either present naturally or systematically produced using the proton or 
electron irradiation. The inhomogeneities both on the microscopic and the mesoscopic scale greatly affect thermody- 
namic and especially dynamic properties of type II superconductors in magnetic field. The magnetic field penetrates 
the sample in a form of Abrikosov vortices, which can be pinned by disorder. In addition thermal fluctuations also 
greatly influence the vortex matter in high T c superconductors, for example in some cases thermal fluctuations will 
effectively reduce the effects of disorder. As a result the H — T phase diagram of the high T c superconductors is 
very complex due to the competition between thermal fluctuations and disorder, and it is still far from being reliably 
determined, even in the best studied superconductor, the optimally doped YBCO superconductor. 1 Difficulties are 
7-H ■ both experimental and theoretical. Experimentally various phases with various (frequently overlapping) names like 
liquid 2 (sometimes differentiated into liquid I and liquid II 3 ), vortex solid, Bragg glass 4 (=pinned solid), vortex glass 
(=pinned liquid=entangled solid, 5 the vortex slush 6 ), were described. To differentiate various phases one should 
understand the nature of the phase transitions between them. Although over the years the picture has evolved with 
various critical and tricritical points appeared and disappeared, several facts become increasingly clear. 

1. The first order 7,8 melting line seems to merge with the "second magnetization peak" line forming the univer- 
sal order - disorder phase transition line. 9,10 At the low temperatures the location of this line strongly depends on 

. disorder and generally exhibits a positive slope (termed also the "inverse" melting 11 ), while in the "melting" sec- 
tion it is dominated by thermal fluctuations and has a large negative slope. The resulting maximum at which the 
magnetization and the entropy jump vanish was interpreted either as a tricritical point 3,12 or a Kauzmann point. 13 
O ■ This universal "order - disorder" transition line (ODT), which appeared first in the strongly layered superconductors 
(BSCCO 9 ) was extended to the moderately anisotropic superconductors (LaSCCO 10 ) and to the more isotropic ones 
like YBCO. 13 ' 14 The symmetry characterization of the transition is clear: spontaneous breaking of the translation 
and rotation symmetry. 

2. The universal "order - disorder" line is different from the "irreversibility line" or the "glass" transition (GT) 
line, which is a continuous transition. 15,16 The almost vertical glass line clearly represents effects of disorder although 
the thermal fluctuations affect the location of the transition. Experiments in BSCCO 17 indicate that the line crosses 
the ODT line right at its maximum, continues deep into the ordered (Bragg) phase. This proximity of the glass line 
to the Kauzmann point is reasonable since both signal the region of close competition of the disorder and the thermal 
fluctuations effects. In more isotropic materials the data are more confusing. In LaSCCO 18 the GT line is closer 
to the "melting" section of the ODT line still crossing it. In YBCO we are not aware of a claim that the GT line 
continuous into the ordered phase. Most of experiments 12 indicate that the GT line terminates at the "tricritical 
point" in the vicinity of the maximum of the ODT line. It is more difficult to characterize the nature of the GT 
transition as a " symmetry breaking" . The common wisdom is that " replica" symmetry is broken in the glass (either 
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via "steps" or via "hierarchical" continuous process) as in the most of the spin glasses theories. 19 

Theoretically the problem of the vortex matter subject to thermal fluctuations or disorder has a long history. An 
obvious candidate to model the disorder is the Ginzburg - Landau model in which coefficients have random components. 
However this model is too complicated and simplifications are required. The original idea of the vortex glass and the 
continuous glass transition exhibiting the glass scaling of conductivity diverging in the glass phase appeared early 
in the framework of the frustrated XY model (the gauge glass). 20,21 In this approach one fixes the amplitude of 
the order parameter retaining the magnetic field with random component added to the vector potential. It was 
studied by the RG and the variational methods and has been extensively simulated numerically. 6 ' 22 In analogy to the 
theory of spin glass the replica symmetry is broken when crossing the GT line. The model ran into several problems 
(see Giamarchi and Bhattacharya in Ref. 23 for a review): for finite penetration depth A it has no transition 24 and 
there was a difficulty to explain sharp Bragg peaks observed in the experiments at low magnetic fields. To address 
the last problem another simplified model had been proven to be more convenient: the elastic medium approach 
to a collection of interacting line-like objects subject to both the pinning potential and the thermal bath Langevin 
force. 25 ' 26 The resulting theory was treated again using the gaussian approximation 4 ' 27 and RG. 21 The result was 
that in 2 < D < 4 there is a transition to a glassy phase in which the replica symmetry is broken following the 
"hierarchical pattern" (in D = 2 the breaking is "one step"). The problem of the very fast destruction of the vortex 
lattice by disorder was solved with the vortex matter being in the replica symmetry broken (RSB) phase and it was 
termed "Bragg glass" . 4 It is possible to address the problem of mesoscopic fluctuation using an approach in which 
one directly simulates the interacting line-like objects subject to both the pinning potential and the thermal bath 
Langevin force. 25 ' 26 In this context the generalized replicated density functional theory 28 was also applied resulting 
in one step RSB solution. Although the above approximations to the disordered GL theory are very useful in more 
"fluctuating" superconductors like BSCCO, a problem arises with their application to YBCO at temperature close 
T c (where most of the experiments mentioned above are done): vortices are far from being line- like and even their 
cores significantly overlap. As a consequence the behavior of the dense vortex matter is expected to be different from 
that of a system of pointlike vortices and of the XY model although the elastic medium approximation might still be 
meaningful. 29 

To describe the non-pointlikc vortices, one has to return to the GL model and make a different simplification. 
One of the most developed schemes is the lowest Landau level (LLL) approximation valid close to the H C 2(T) line. 30 
Such an attempt was made by Dorsey, Fisher and Huang 31 in the liquid phase using the dynamic approach 32 and 
by Tesanovic and Herbut for columnar defects in layered materials using supersymmetry. 33 It is the purpose of this 
paper to study the glass transition using the replica formalism. We quantitatively study the glass transition in the 
same model, with the disorder represented by the random component of the coefficients of the GL free energy. The 
most general hierarchical homogeneous (liquid) Ansatz 34 and its stability are considered to obtain the glass transition 
line and to determine the nature of the transition for various values of the disorder strength of the GL coefficients. 
Then we place the glass line on the phase diagram of YBCO and compare with experiments and other theories. 

The paper is organized as follows. The general disordered GL model is introduced in section II and the gaussian 
variational replica method is presented in section III. Next we study in some detail the model either with \tp\ disorder 
in Section IV, with less details the \ip\ disorder in Section V and obtain the phase transition lines in those two cases. 
In Section VI, the general model containing both the |?/>| 2 disorder and the |?/>| 4 disorder is treated briefly. In Section 
VII, we compare our results with the experimental data and conclude in section VIII by summarizing our results. 



II. DISORDER EFFECTS IN THE GINZBURG - LANDAU DESCRIPTION OF THE TYPE II 

SUPERCONDUCTOR 



A. Ginzburg - Landau free energy 

We start from the Gibbs energy of the ideal homogeneous sample (no disorder): 



G 



2m* ± 



Dtp 



(H-Bf 



+ a » + -(v>*V) 2 + . 



(1) 



Here a' = a(T — T c ) and b' are constant parameters, D = (—i~h\7 + ^- A) is the covariant derivative, A is the vector 
potential, the magnetic field B = V x A , H is the external magnetic field, and are the effective masses in 
directions perpendicular and parallel to the field respectively. Mesoscopic thermal fluctuations are accounted for via 



3 



Boltzmann weights 



cxp <^ = > (2) 



The model provides a good description of thermal fluctuations as long as 1 — t — b « 1, where t = b = Bj H c2 

(h = H/H C 2), H C 2 — <I>o/(27r£ 2 ) and £ is the coherence length. In this case the higher order terms like |?/>| 6 can be 
omitted (detail notation can be found in Ref. 35). The 3D GL model describes materials with not too high anisotropy 
(for a recent evidence of validity of this assumption in YBCO see Ref. 36). In strongly anisotropic materials, a 
model of the Lawrence - Doniach type is more appropriate. 37 

Within the GL approach the pointlike quenched disorder on the mesoscopic scale is described by making each of 
the coefficients of the mesoscopic GL energy a random value centered around a certain constant value given in eq. 
(2). For example effective masses can be disordered 



m*- 1 -» m*- 1 (1 + U(x)) ; (3) 



U(x)U(y) = PS(x-y). 

The parallel effective mass mfj might also have the random component which we neglect (it is relatively small since 
to*| is typically very large), though it can be incorporated with no additional difficulties. This type of disorder is 
sometimes called the SI disorder since it originates in part from the inhomogeneity of the electron mean free path 
I in Gor'kov's derivation. From the BCS theory, effective mass is m* = 2m e ^1 + 16 g^^y ; ^ in the clean limit and 

m* = 2to c 7 -^$rf m * ne dirty limit. Relation to the notations of Ref. 37 (chapter II in this reference) is the following: 
U is —Sm a b/m a b and P — 7 m /TO 2 ib . Note however that, in addition to the random distribution of I, disorder in vp 
and T c (the density of states and interaction strength) can also affect to*. 

The other two parameters in the GL equations are a = an d b = 7 ^^ <iF (tf) ■ The coefficient of the 

quadratic term is called ST disorder since it describes a local deviation of the critical temperature. Introducing a 
random component in |?/>| 2 term: 



a' -» a'{\ + W(x)); W(x)W(y) = RS(x - y). (4) 

In notations of Ref. 37 the random field W{x) = —S a /a, R = j a /a 2 . When thermal fluctuations of the vortex degrees 
of freedom can be neglected, these two random fields would be sufficient (they control the two relevant scales £ and 
A). The reason is that one can set the coefficient of the third term \ip\ to a constant by rescaling. However in the 
presence of thermal fluctuations the coefficient of \tp\ 4 also should be considered as having a random component. It 
cannot be "rescaled out" since it affects the Boltzmann weights. We will see later that at least within the lowest 
Landau level approximation this term is crucial in inducing certain glassy properties of the vortex matter state. We 
therefore introduce its disorder via 



b'^b'(l + V(x))- V(x)V(y)=QS(x-y). (5) 

In unconventional superconductors, even without disorder, the phenomenological GL model has not been reliably 
derived microscopically. The coefficients and their inhomogencities therefore should be considered as phenomeno- 
logical parameters to be fitted to experiments. We assume that U, R and Q have weak dependencies on field and 
temperature. The assumption of the weak temperature and field dependence of the disorder strengths U, R and Q, 
as that of any parameter in the GL approach, should be derived in principal from a microscopic theory assuming 
random chemical potential or should be justified by fitting to experiments. For simplicity the white noise distribution 
is considered 



p[u, n .r cxp ^ 



L 



U{x) 2 t W{x) 2 V{xf 



2P 2R 2Q 

for random components. The free energy of superconductor after averaging over the disorder is 

T 



F = — 



norm 



f p[U,W,V]log 
Ju.w.y 



1 



exv[-9bP}-fdis[U,W,V,: 
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(6) 
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/ 7^ U(x) 

L 2ml 





.9 = G/T- f dis [U,W,V^} = ^ / £^U(x) DV +a'W(xM 2 + -V(xM 4 , (7) 
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where norm = J uwv p[U, W, V] is a normalization factor. 

To make the physical picture clear, we rescale the coordinates as x — > £x, y — > £y, z — > — with anisotropy parameter 

1 /2 

defined by 7 = ^m*|/m^^ . The order parameter is scaled as i[> 2 — > 2 jp L tp 2 . The dimensionless free energy becomes 
simpler looking: 



t — 1 . , l2 1, ,,4 K 2 (b-hf 
+ ^-|^| 2 + ^l^| 4 + V 4 



(8) 



where cj = \/2Gm 2 t. The Ginzburg number is Gi = 32 [7rA 2 T c 7/ ($q0] > where A is the magnetic penetration depth. 
The last term can be ignored in calculating F as the k is very big in high T c superconductors and the last term is 
order of Similarly the random component and the distribution become: 



f dis [u,w,v^\ = i j{-hj{x)r&il>- 



t —}-W{x)\^\ 2 



-V{x)\ 



p[U,W,V] = exp 



U(x) 2 W{x) 2 V{xf 



(9) 
(10) 



7 j x V 2P 2R 2Q 

The model however is highly nontrivial even without disorder, and to make progress, further approximation is needed. 



B. Lowest Landau level approximation 



The lowest Landau level (LLL) approximation 30 is based on constraint —D 2 tjj = bip. Over the years this model 
has been studied by various methods, analytic and numerical. 35 ' 38,39 The (effective) LLL model is applicable in a 
surprisingly wide range of fields and temperatures determined by the condition that the relevant excitation energy e 
is much smaller than the gap between Landau levels 2heB/(cm±). 13 

/ 5 / 2 \ 1 / 3 / 5 / 2 \ 2 / 3 

The free energy after further rescaling x — > x/V~b,y — > y/Vb,z — > z ( 2 -^j L ) , ip ~ > ( 2 fea , 7r ) simplifies 

within the LLL approximation to: 



Jlll 



2 5 / 2 



7T 



d 3 a 



(11) 



Not surprisingly the number of independent constants in LLL is one less than in the general model. This fact leads 
to the "LLL scaling" relations (of course the disorder terms will break LLL scaling). As a result the simplified model 
without disorder has just one parameter - the (dimensionless) scaled temperature: 



/ 9 \ 2/3 



The disorder term becomes: 



2 + -V(x 



in which only combination of W and U enters Q(x) = | 
gaussian 



2(t-l)(t) 2/3 W(x)-bU(x) 



p(Q, V) = exp 



n{ x ) 2 v{xf 



r n(x) 

L 2r' 



2q 



with two variances 



r = 



Q = 



1/3 
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V2e 



1/3 



b 2 p 



(12a) 
(13) 

. Its distribution is still 
(14) 

(15) 



2^ Q - 
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To treat both the thermal fluctuations and disorder we will use the replica method to integrate over impurity distri- 
bution followed by gaussian approximation. 



III. REPLICA TRICK AND GAUSSIAN APPROXIMATION 

A. Replica trick 

We will use the replica trick to evaluate the disorder averages. The replica method is widely used to study disordered 
electrons in the theory of spin glasses, 19 disordered metals and was applied to vortex matter in the London limit. 27 ' 40 
Applying a simple mathematical identity to the disorder average of the free energy one obtains: 



F = -T\im—(Z n - 1). (16) 

n^On 

The averages of Z n is the statistical sum over n identical "replica" fields ip a , a — 1, ...,n: 

W=—[ p[Q,V\T[f exp{~f[^ a }-f dls [n,V^ a }}. (17) 
norm J ny a J-4> a 

The integral over the disorder potential is gaussian and results in: 



Z n = / exp 



b 



(18) 



fab = r'\^ a \ 2 \^ b \ 2 + q -{^a) 2 {rM 2 - 

This model is a type of scalar field theory and the simplest nonperturbative scheme commonly used to treat such a 
model is gaussian approximation. Its validity and precision can be checked only by calculating corrections. 

B. Gaussian approximation 



We have assumed that the order parameter is constrained to the LLL and therefore can be expanded in a basis of 
the standard LLL eigenfunctions in Landau gauge: 

V>» = norm j e ^ +x ^ exp {-Uy + k) 2 \^ a {k). (19) 



We now apply the gaussian approximation which has been used in disorder in the elastic medium approach, 27 ' 40 
following its use in polymer physics. 41 The gaussian approximation was applied to the vortex liquid within the GL 
approach in 30,39 The gaussian effective free energy is expressed via variational parameter 35,41 fi ab which in the 
present case is a matrix in the replica space. The correlator is parametrized as 

(r a (k, k z )M~k, -k z )) = G ab (k z ) = k2 25/2?r (20) 

The bubble integral appearing in the free energy is very simple: 

V^2 f 1 

(tp* a {x, y, z)TP b (x, y, z)) = — / — — = 2^ = 2m ab . 
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As a result the gaussian effective free energy can be written in a form: 



n feff = E||^/ [LogG-\k z )+(^+a T ^G{k z )-I 

~ E | 23~/2^ r ' l mQf »| 2 + ~<f (l m ab| 4 + ^m aa m bb \m ab \ 2 ^j 

a,b K 

= 2^ |m qq + a T m aa + 2 (m aa ) 2 | 

a 

-2 ^ |r |m afc | 2 + q Q |m afc | 4 + m aa m bb \m a . 



+ 4(m aQ ) 2 



(21) 



I 2 

where we discarded an (ultraviolet divergent) constant and renormalization of ar and rescaled the disorder strength: 
r = W ^r, q=—q- 

We start with a simple case in which only the \ip\ type of disorder is present. More precisely we take q = 
and return to the general case in section IV. This model has been already discussed using different method (the 
Sompolinsky dynamic approach) in the unpinned phase in Ref. 31. 

IV. NONZERO EDWARDS - ANDERSON ORDER PARAMETER AND ABSENCE OF THE REPLICA 
SYMMETRY BREAKING WHEN ONLY THE \tp\ 2 DISORDER IS PRESENT 

A. Hierarchical matrices and impossibility of the continuous replica symmetry breaking 

In this section we neglect the \ip\ 4 disorder term. It is convenient to introduce real (not necessarily symmetric) 
matrix Q abl which is in one to one linear correspondence with Hermitian (generally complex) matrix m ab via 

Q ab = rc[m ab ] + im[m a6 ]. (22) 

Unlike m abl all the matrix elements of Q ab arc independent. In terms of this matrix the free energy can be written as 



\USS = E { K^aa + ^Qaa + 2 (Q aa f} ~ r £ Q* 6 . (23) 
a a.b 

Taking derivative with respect to Q ab gives the saddle point equation for this matrix element: 

J-Jt- = -\ [(1 - ») ( m ~ 2 ) ab + CC] + a T 8ab + ^QaaSab - 2rQ ab = 0. (24) 

Since the electric charge (or the superconducting phase) U{1) symmetry is assumed, we consider only solutions with 
real m ab . In this case m ab = Q ab is a symmetric real matrix. General hierarchical matrices m are parametrized using 
the diagonal elements fh and the Parisi's (monotonically increasing) function m x specifying the off diagonal elements 
with < x < l. 41 Physically different x represent time scales in the glass phase. In particular the Edwards - Anderson 
(EA) order parameter is m x= i = M > 0. 

A nonzero value for this order parameter signals that the annealed and the quenched averages are different. The 
dynamic properties of such phase are generally quite different from those of the nonglassy M = phase. In particular 
it is expected to exhibit infinite conductivity. 20 ' 31 We will refer to this phase as the "ergodic pinned liquid" (EPL) 
distinguished from the "nonergodic pinned liquid" (NPL) in which, in addition, the ergodicity is broken. 

However in the present model RSB does not occur. In terms of Parisi parameter fh and m x the matrix equation 
eq.(24) takes a form: 

-nr^ 2 + a T + (4 - 2r) fh = (25) 



(m 2 ) x + 2rm x = 0. 
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Dynamically if m x is a constant, pinning does not results in the multitude of time scales. Certain time scale sensitive 
phenomena like various memory effects 42 and the responses to "shaking" 11 are expected to be different from the case 
when m x takes multiple values. If m x takes a finite different number of n values, we call n — 1 step RSB. On the 
other hand, if m x is continuous, the continuous replica symmetry breaking (RSB) occurs. 

In order to show that m x is a constant, it is convenient to rewrite the second equation via the matrix fj,, the matrix 
inverse to m: 

(H 2 ) x + 2r(p- 1 ) x = 0. (26) 
Differentiating this equation with respect to x one obtains; 



M x -r(M x r 2 ]^=0, (27) 
where we used a set of standard notations in the spin glass theory: 41 

W^M-W-Nii (/*x>=/ dx V x ; [/*]*=/ dy(n x -v y ). (28) 

Jo Jo 

If one is interested in a continuous monotonic part ^ 0, the only solution of cq.(27) is 

M^r 1 ' 3 (29) 

Differentiating this again and dropping the nonzero derivative -4^ again, one further gets a contradiction: = . 
This proves that there are no such monotonically increasing continuous segments. One can therefore generally have 
cither the replica symmetric solutions, namely m x = M or look for a several steplike RSB solutions. 19 We can show 
that the constant m x solution is stable. Therefore, if a stcplikc RSB solution exists, it might be only an additional 
local minimum. We explicitly looked for a one step solution and found that there is none. 

B. Two replica symmetric solutions and the third order transition between them 

1. The unpinned liquid and the "ergodic glass" replica symmetric solutions 

Restricting to RS solutions, m x = M, the saddle point equations eq.(25) simplify: 

-e~ 2 + (or + 4m) - 2re = 0; (30) 
M (e~ 3 - r) =0, 



where e = m — M. Energy of such a solution is given by 

^ e// - 26" 1 - 2 - £ - 2 M + 2a T m + Am 2 - 2r (e 2 + 2eM) . (31) 



The second equation eq.(30) has a replica index independent (diagonal) solution M = 0. In addition there is a non 
diagonal one. It turns out that there is a third order transition between them. 
For the diagonal solution e = rh and the first equation is just a cubic equation: 

— in 2 + (a T + 4m) - 2rm = 0. (32) 
For the non diagonal solution the second equation gives e = r 1//3 , which, when plugged into the first equation, gives: 

m = \ (3r 2 / 3 - or) ; M - \ (ir 2 '* - a T ) - r^ 3 . (33) 

The matrix m therefore is 

m ab = r-^Sab + M, (34) 
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which results in the following value of the free energy: / = Qr 1 ^ 3 — j (3r 2//3 — ay) 2 . 
The two solutions coincide for 

a T = r~ 1/3 (3r - 4) . (35) 

Since in addition to the energy, the first and second derivatives of the energy, — 2r -1 / 3 and ^Jr — — \ respectively, 
coincide (the fourth derivatives are different though), the transition is a third order one. 

2. Stability domains of the two solutions 

In order to prove that a solution is stable beyond the set of replica symmetric matrices m one has to calculate the 
second derivative of free energy (called Hessian in Refs. 19,43) with respect to arbitrary real matrix Q ab defined in 
eq.(22): 

;2 



n S f eff 
2 
1 



H(ab)(cd) = TT TTi xn ( 36 ) 

2 oQabdQcd 



^[K 2 ) QC K 1 ) db -K 2 ) a ,(™- 1 )J + 

+4:6 ac 5 bd S ab - 2rd ac S bd . 

We will use a simplified notation for the product of the Kronecker delta functions with more than two indices: 
Sac^bdSab = Sabcd- For the diagonal solution the Hessian is a very simple operator on the space of real symmetric 
matrices: 

H(ab)(cd) — Cllabcd + CjJabcd, (37) 

where the operators / (the identity in this space) and J are defined as 

I = Sachd', J = Sabcd (38) 

and their coefficients in the diagonal phase are: 



Cl = 2 (m~ 3 - r) , cj = 4 (39) 

with to being a solution of cq.(32). The corresponding eigenvectors in the space of symmetric matrices are V( cd ^ = 
AS c d + B. To find eigenvalues A of if we apply the Hessian on V. The result is (dropping terms vanishing in the limit 
n->0): 

H(ab){cd)Vcd = A (a + cj) S ab + B (a + cjS ab ) = X (AS ab + B) (40) 

There two eigenvalues: = cj and A^ 2 ^ = cj + cj. Since cj = 4 > 0, the sufficient condition for stability is: 

a = 2 (to~ 3 - r) > 0. (41) 

It is satisfied everywhere below the transition line of eq.(35), see Fig l.( ar, r phase diagram). The analysis of stability 
of the non diagonal solution is slightly more complicated. The Hessian for the non diagonal solution is: 



where new operators are 



H(ab)(cd) = c-vV + CuU + CjJ, 

V(ab)(cd) = S ac + 5bd \ U( a b)(cd) = 1 



(42) 
(43) 
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and coefficients are 

c v = -3Mr 2/3 ; cjj = AM 2 r 1/3 ;cj = 4 (44) 

In the present case, one obtains three different eigenvalues, 19 ' 43 A^ 1 ' 2 ^ = 2 ^1 ± \f\ — 4Mr 2 / 3 j and A^ = 0. Note 

that the eigenvalue of Hessian on the antisymmetric matrices are degenerate with eigenvalue A^ in this case (we 
will come back later on this eigenvalue). For M < the solution is unstable due to negative X^ 2 \ For M > 0, both 
eigenvalues are positive and the solution is stable. The line M — coincides with the third order transition line, 
hence the non diagonal solution is stable when the diagonal is unstable and vise versa. We conclude that there is no 
glass state in the vortex liquid without the \ip\ 4 disorder term. The transition does not correspond to RSB. Despite 
this in the phase with nonzero EA (NEA) M order parameter there are Goldstone bosons corresponding to A*- 3 -* in 
the replica limit of n — > 0. The criticality and the zero modes due to disorder (pinning) in this phase might lead to 
great variety of interesting phenomena in statics and dynamics. These have not been explored yet. However, as we 
show in the next section, the random component of the quartic term changes the character of the transition line: the 
replica symmetry is broken on the one side of the line. For simplicity in the next section, we consider first a case with 
a random component of \ip\ and no random component of \ijj\ , and return to the general case in section V. 

V. THE GLASS TRANSITION FOR THE \ib\ 4 DISORDER 



A. Continuous replica symmetry breaking solutions 

In this section we neglect the r \ip\ 2 term disorder. Although it is always present, as we have seen in the previous 
section, at least within the gaussian approximation, it docs not cause replica symmetry breaking. Therefore one 
expects that although it certainly influences properties of the vortex matter, for example, the melting transition line 
to lower fields and temperatures, 13 it's role in qualitative understanding of RSB effects is minor. The only other 
disordered term within the LLL approximation considered in this paper is the \ip\ 4 disorder term. As was discussed 
in section II, at least within the BCS theory, it is expected to be smaller than the \tp\ 2 disorder, g«r. Even it could 
be very small, however, as we show here, it leads to qualitatively new phenomena in vortex matter. The r = free 
energy after integration over k z becomes: 



feff = E {( m X )aa + a T m aa + 2 {m aa f} - \ \ \™>abf + m aa m bb |m afc | 2 j 



(45) 



a,b 

In terms of the real matrix Q a ^ denned in eq.(22) the free energy can be written as 



\Utt - J2{( m ^)aa+ a TQaa + 2(Q aa ) 2 } (46) 
~ a 

-1 E (Wab + \QlbQla + QaaQbbQlb) (47) 



a,b 



Taking a derivative with respect to Q a b gives the saddle point equation for this matrix: 
n Sf 



2 SQ ab 



~\ ~ ( m 2 )q6 + CC ] + a T S ab + 4:Qaa8ab 
-1 {\Qlb + \QabQl a + IQabQaaQbb + S ab £ e Q ee (Q 2 £ + Ql a )) 



= (48) 



Using the hierarchical symmetric matrix parametrization of its symmetric part (the antisymmetric will not be impor- 
tant for most of our purposes), it takes a form 

-m^ 2 + a T + 4m-q (Zm 3 + 2m 2 fhj = (49) 
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(m- 2 ) x + q(m 3 x + 2m 2 m x ) = 0. (50) 
As in the previous section, it is convenient to rewrite the second equation in terms of fi, the inverse matrix of m: 



Differentiation of this equation with respect to x leads to: 



0. 



(51) 



{2 {fJ,} x - q [3 ((p-%) 2 + 2m 2 ] ({ M } J" 2 } x^ = 0. 



(52) 



For a continuous segment 7^ one solves eq.(52) for (/j, 1 ) x in terms of {n} x getting now a more complicated 
result: 



1 '(WJ -™ ! 



(53) 



Differentiating this equation with respect to x again one obtains: 



(54) 



Instead of solving this for {^} x we present x as function of {n} x ~- 



x = 



3q 



2(M,f 



(55) 



Thus the solution will be given by eq.(55) in the segment if 7^ and constant fj, x in the other segments. In 
principle this would allow for a numerical solution. One could actually solve the equation near the transition line 
using the method in Ref. 34. The situation is completely different compared to that of the \ip\ 2 disorder. In the 
present case a stable RSB solution exists. We will turn first however to the replica symmetric solutions and determine 
their region of stability. In the unstable region of the replica symmetric solutions, the RSB solution of eq.(55) will be 
the relevant one. 



B. Two replica symmetric solutions 



1. Solutions 



Here we briefly repeat the steps leading to the RS solutions for the disorder omitting details. The saddle point 
equations eq.(25) for the RS matrices m x = M are: 



-e~ 2 + (a T + Am)~q [5m 3 - M 3 - 2M 2 m - 2m 2 M] = 



(56) 



M [2e~ 3 - g (M 2 + 2m 2 )] =0 
where e = fa — M. Energy of such a solution is given by 



2 = e- 1 - e~ 2 M + a T m + 2m 2 - | (5m 2 - M 2 ) (m 2 + M 2 ) 



(57) 



(58) 
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For the diagonal solution M = 0, s = m and the first equation takes a form: 



-to 2 + a T + 4m - 5qm 3 = 0. (59) 

The non diagonal solution in the present case is more complicated, but the condition determining the transition line 
between the two (equivalently the appearance of the nonvanishing EA order parameter) is still very simple: e = q~ x ^ 
as M = on the line. Along the line the scaled temperature is: 

4 = 2(3<? 2/5 -2q- 1 / 5 ). (60) 

It is still the third order transition line similar to the \i/;\ 2 disorder case and one has zero modes in NEA sector, while 
no such modes exist in the M = phase. The stability analysis with respect to configurations which are replica 
symmetric however gives completely different results compared to that of the \if>\ 2 disorder. 

2. Stability region of the RS solutions. 
The Hessian now has several additional terms 

n 6 2 f 

H( - ab ^ = 2 Sm ab 5m cd = ^c( m ~ )db+( m ~ ) ac ( m " ) rfb + 4 <Wrf (61) 

(^S ac SbdQab + 3 {^ac^bdQba + ^ad^bcQabQba) 
2(5 ac SbdQaaQbb + SbcdQabQaa + 5 acdQ abQ bb) 
+ Sab [Scd (Qlc + Qla) + 2 (SacQadQdd + SadQcaQcc) 

For replica symmetric solutions Q a b = m-ah = £d a b + M the Hessian can be represented as 

H = c+I+ + c-l- + CuU + c v V + cjJ + c K K + c N N, (62) 
where new operators I± , K, N are defined as 

I± = ^ {5 ac 5bd ± S ad 6 bc ) ; K = S ab S cd ; N = S abc + S abd + S acd + S bcd . (63) 

The coefficients are 



c+ = 2e- 3 - q (3M 2 + 2m 2 ) ; (64) 

c_ = 2e~ 3 - q (M 2 + 2m 2 ) ; 

c v = 4M 2 £- 5 ;c v = -3M £ - 4 ; 

Cj = 4 - q [5 (to 2 - M 2 ) + 8fh (m - M )] ; 

= -2qM 2 ; c N = -2qrhM. 

Generally the Hessian have four different eigenvalues: 43 



A (1 ' 2) = c+ + - [c,, + Ac N ± Vcj (cj + Sc v + 8 CJV )) ; A (3) - c+, A (4) = c_ (65) 

Note that there are new matrices like I+,I- when q ^ 0. In the case of q = 0, c + = c_, so that only operator 
/ = /++/_ appears in this case. Actually there is which is the eigenvalue of Hessian on the antisymmetric 
matrices. However = c > is always hold on the RS solutions so that it can be ignored in determining the 
instabilities of those RS solutions. Since the stability analysis is quite complicated, we divide it into several stages of 
increasing complexity. 



12 



3. Stability of the states on the diagonal - off diagonal "transition" line 

The easiest way to see that the RS solutions can be unstable is to look first at the transition line dj>, eq.(60). On 
the transition line one has 



and the eigenvalues simplify to 
Therefore it is unstable for q > q* 

marginally stable at a single point 



c± = cu = c v = c N = 0; cj = 4 - 13«? 3/5 ; (66) 

A (3) =Q; A (l,2) = 4 _ 13q 3/5 (6?) 



4 \ 



5/3 



4 = -| (f) V3 --3.2 (69) 

and stable for q < q l . We studied numerically the stability on both sides of this line, see Fig. 2. The diagonal (liquid) 
solution is stable below the line (a^ = 2(3q 2 ^ 5 — 2q~ 1 / 5 )) for q > q t . The line when q > q l the phase transition line 
(liquid to glass) is changed to a different line which will be discussed in the next subsection. 

4- Stability of the diagonal solution 

Equation for rh, coefficients in Hessian and eigenvalues are: 

-m~ 2 + a T + 4m - 5qm 3 = 0; (70) 

c + = 2to~ 3 -2qm 2 ; Cj = 4-I3qm 2 ; (71) 
A (1 ' 3) = c+; A (2) = c+ + cj = 4 + 2m" 3 - 15qm 2 . 

While A*- 1 ' is positive, A^ 2 -* is positive only below the line defined parametrically via 

5 - 8m 3 4m 3 + 2 

aT = ^ 2 ~' ^^5^' (72) 

marked by dotted line in Fig. 2. and this line is the phase transition line (liquid to glass) when q > q l . The former 
diagonal - off diagonal line above the tricritical point is not a phase transition and is left as a light dashed line to 
show that slope of the line below tricritical point and that of the real transition line is different, ft turns out that the 
line of eq.(72) is a transition line into a RSB state, namely the irreversibility line. 

5. Stability of the off diagonal solution 

The equations take a form: 

-a' 2 + a T + 4m - q [5m 3 - M 3 - 2M 2 m - 2m 2 M] = (73) 

2e~ 3 - q (M 2 + 2m 2 ) = 0. (74) 

The coefficients in the expansion of Hessian are: Unlike the case of the IV^ 2 disorder, A^ 1 ^ < for each such a solution. 
Therefore the diagonal state directly goes over into a RSB glass state. It follows however two lines. The eq.(72) above 
the tricritical point and eq.(60) below it, see Fig. 2. 



VI. RSB IN THE GENERAL CASE 
A. General hierarchical gaussian variational Ansatz 



The free energy 
/ 



= {( m l ) aa + a rm aa + 2(m aa ) 2 } - ^ \2rm 2 ab - |(to* 6 + 4m aa m bb m 2 a 



2n 

a a.b 

lead on the replica symmetric sector to the following equations: 

-e~ 2 + a T + 4m - 2re - q [5m 3 — M 3 — 2M 2 m - 2m 2 M] = 0, 

M [2e~ 3 -2r-q (M 2 + 2m 2 )] = 0. 

For the diagonal solution M — one obtains 

—wT 2 + ax + 4m — 2rm — 5qm 3 = 0. 

The off diagonal solution on the bifurcation line obeys 

m~ 3 — r — (jrn 2 = 0. 

Hessian for the general RS solution takes a form of eq.(62) with coefficients 

c+ = 2e~ 3 — 2r — q (3M 2 + 2m 2 ) ; c_ = 2e~ 3 - 2r - q (M 2 + 2m 2 ) ; 

Cu = 4M 2 e~ 5 ; c v = -3Me" 4 ; c K = -2qM 2 ; 

c N = -2qrhM; cj = 4 - q [5 (m 2 - M 2 ) + 8m (m - M)] . 

On the bifurcation line it simplifies: 



c± = 2m 3 — 2r — 2qm 2 ; Cu = cy = Ck = Cn = 0; cj = 4 — 13qm 2 . 



The eigenvalues are 



A^=c + ; A( 3 >= C++CJ . 
Therefore the Hessian vanishes A^ 1 ' 2 ^ = A^ = for the tricritical (branch) point defined by 

Stability of the diagonal solution Hessian and eigenvalues in this case are: 

c ± = 2m" 3 -2r- 2qm 2 ; cj = A- 13<?m 2 ; 
A (1 ' 2) = c±; A (3) =c ± +c./ = 4 + 2m" 3 -2r-15 g m 2 . 

Below the tricritical point we solve equation = perturbatively in q : 

m = r- 1 / 3 (l-|r 5 / 3 ) + 0(q 2 ) 

and substitute m into eq.(77) to determine the "weak disorder" part of the glass transition line: 

<7i 4 - r f 5 47-4/3 \ ^ 2^ 



V3 ^ r 3 y 
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Above the tricritical point, namely for larger disorder, one solves the equation A*- 3 -* = perturbatively in q around 
the tricritical point of eq.(82), q l {r) = ^ (r + j^) 3 ^ 2 '■ 

• a =( r+ s)" /3 ( 1 -24Ti3; a ) +0 < a! > ; a = ?f?)- L (86 » 

The NEA RS solution is unstable everywhere as c + < (c + < c_ = 0). We therefore obtain the glass transition 
line with RSB in the general case. To compare it with experiment one has to specify phenomenologically the precise 
dependence of the GL model parameters on temperature. Next section is devoted to this. 

VII. COMPARISON OF THE THEORETICAL AND THE EXPERIMENTAL PHASE DIAGRAMS 

A. General picture and comparison with other theories 



As was discussed in Introduction, the interplay between disorder and thermal fluctuations makes the phase diagram 
of high T c superconductors very complicated. As a result of the present investigation together with the preceding 
one 13 regarding the order - disorder transition the following qualitative picture of the T — H phase diagram of a 
3D superconductor at temperatures not very far from H C 2(T) (so that the LLL approximation is valid) emerges, see 
schematic diagram Fig. 3. There are two independent transition lines. 

1. The positional order- disorder line. 

The unified first order universal order - disorder line comprising the melting and the second peak sections separates 
the homogeneous and the crystalline (Bragg) phases. The transition is therefore defined by the translation and 
rotation symmetry only, and the intensity of the first Bragg peak can be taken as order parameter. The broken 
symmetry is not directly related to pinning, however the location of the line is sensitive to disorder. One sees on 
Fig. 3 that as the disorder strength n increases the melting line (solid line) curves down at lower point merging with 
the second peak segment. The effect of disorder is quite minor in the high temperature region, in which the thermal 
fluctuations dominate, but become dominant at low temperatures. The line makes a wiggle near the experimentally 
claimed critical point. The "critical point" is reinterpreted 13 as a (noncritical) Kauzmann point in which the latent 
heat vanishes and the line is parallel to the T axis in the low temperature region. The surprising "wiggle", which 
appears in 3D only, has actually been observed in some experiments. 14 It is located precisely in the region in which the 
thermal fluctuations and the disorder compete. One might expect that the line just has a maximum like in BSCCO, 
but the situation might be more complex. Thermal fluctuations, on the one hand side, make the disorder less effective, 
the less disorder effect will favor solid, but, on the other hand side, they themselves melt the solid. The theoretical 
magnetization, the entropy and the specific heat discontinuities at melting line 13 compare well with experiments 3,12 . 
The low temperature segment of the disorder dominated second peak line is weakly temperature dependent. Its field 
strongly decreases as the disorder strength increases. 

2. The glass transition line. 

The second line is the glass transition line, the dashed line in Fig. 3. We assume very small q and draw on Fig. 3 
the glass lines for three different values of n. One observes as expected, that, as the disorder strength increases, the 
line moves towards higher temperatures. We have not calculated the glass line in the crystalline state, but anticipate 
that it depends little on the crystalline order. This is consistent with observations made in Ref. 33 in which it was 
noticed (in a bit different context of layered materials and columnar defects) that lateral modulation makes a very 
small difference to the glass line although it is obviously very important for the location of the order - disorder line. 
So we just continue the glass line in the homogeneous phase into the crystalline side (still marking it by a dotted 
lines). If the glass lines of the liquid side and solid side join to a single glass line, then the glass line must cross the 
ODT line right at the Kauzmann point. The theoretical calculated intersection point in Fig. 3 marked by black blob 
do not appear at the Kauzmann point. We attribute this to different approximations made to draw the ODT and 
GT line. That the glass line crosses the ODT line right at the Kauzmann point has been indeed supported by some 
evidences in an experiment in BSCCO 11 , thought two lines crossing exactly at the Kauzmann point or not is still an 
open question experimentally and theoretically 

Consequently there are four different phases, see Fig. 4: pinned solid (= Bragg glass), pinned liquid (=vortex 
glass) , unpinned solid (or simply solid) and unpinned liquid (or simply liquid) . The four phases of the vortex matter 
are also expected to be present in the case of layered quasi 2D superconductors and it was shown by in BSCCO in 
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17 that the conductivity and the magnetization both indicate the glass line crossing the ODT line near the maximum 
(Kauzmann point). There is evidence of the crossing of the ODT and GT line also in LaSCCO 18 and YBCO 12 ' 16 , 
but the line seems to lie very close to the melting line (the segment of ODT line in the high temperature). In Ref. 
16 a great care was paid to distinguish the two lines, and it was found that at low field the GT line crosses again the 
melting line so that the lower field segment of the melting line is again separating two glassy or pinned states. 

Now we compare the results on the phase diagram with other theories. The only known theory providing both the 
ODT and the glass lines is the density functional model. The picture advocated in Ref. 28 on basis of the replica 
density functional calculation in the framework of thermodynamics of pinned line objects is qualitatively different from 
the present one. In this theory the glass line does not intersect the ODT line and therefore there is no unpinned solid 
phase. The comparison of the theory studied in this paper with the replica density functional theory is complicated 
by the fact that the applicability ranges of two theories are different. In addition there are other phenomenological 
approaches, most notably that of ref. 44 , based on Lindermann criterion to map the ODT line. 



B. Application to optimally doped YBCO 



We plot the ODT and the glass transition line in Fig. 5 along with the experimental data on optimally doped 
YBCO of Ref. 12. The theoretical lines use the fitting parameters of ref. 13 which fitted the ODT line and 
assume the \ip\ disorder strength q is small. We believe the value of this parameter should be measured directly 
using a replica symmetry breaking dynamical phenomenon rather than fitted in thermodynamics, though it might be 
possible to adjust q, so that the curve can fit the glass transition line better than the case of q = 0. As was argued 
in 13 (and found to be consistent with experiments on the ODT line and the magnetization, entropy and specific 
heat discontinuities at melting line) that the general dependence of the disorder strength on temperature near T c is: 
r(t) = n(l — t) 2 /t with n = 0.12. The formula interpolates the one used at lower temperatures in ref. 37 with our 
dimcnsionlcss pinning parameter n (proportional to the pinning centers density) related to the "pinning strength" 
of Blatter et al 37 by n = 77^/7rG« 1 / 2 £ 3 . Note that experimentally the glass transition line at lower fields is not 
measured precisely. Different experiments locate it at various places using different criteria. As mentioned above 
experimentally the GT line often crosses the melting line again at very low fields. This question perhaps cannot 
be addressed within the lowest Landau level approximation valid at fields above IkG. The ODT line in the melting 
region is very well established experimentally by great variety of techniques. However the "second peak" segment is 
only poorly determined due to difficulty to define it in the essentially dynamic magnetization loops analysis. Recently 
developed muon spin rotation and neutron scattering methods might be very helpful in that respect. 



VIII. SUMMARY 



To summarize we considered the effects of both thermal fluctuations and disorder in the framework of the GL 
approach using the replica formalism. Flux line lattice in type II superconductors undergoes a transition into three 
"disordered" phases: vortex liquid (not pinned), homogeneous vortex glass (the pinned liquid or the vortex glass) and 
the Bragg glass (pinned solid) due to both thermal fluctuations and random quenched disorder. We show that the 
disordered Ginzburg - Landau model (valid not very far from H C 2) in which only the coefficient of a term quadratic 
in order parameter ip is random, first considered by Dorsey, Fisher and Huang, leads to a state with nonzero Edwards 
- Anderson order parameter, but this state is still replica symmetric. Namely there is no ergodicity breaking and 
no multiple time scales in dynamics are expected. However, when the coefficient of the quartic term in ip in the GL 
free energy also has a random component, replica symmetry breaking effects appear (with ergodicity breaking). The 
location of the glass transition line in 3D materials is determined and compared to experiments. The line is clearly 
different from both the melting line and the second peak line describing the translational and rotational symmetry 
breaking at high and low temperatures respectively. The phase diagram is therefore separated by two lines into four 
phases mentioned above. In principle we could obtain the RSB solution near the phase transition line by expanding 
the equations around the phase transition line as in the spin glass theory, see, for example ref. 34, and we found that 
the RSB is continuous. Thus RSB states involve multiple time scales in relaxation phenomena. 

It is natural to expect and is confirmed that the glass (irreversibility) line crosses the "order - disorder" line not 
very far from its Kauzmann point. We are not sure if the crossing is right at Kauzmann point. If two GT lines (on 
liquid side and solid side) are joined, the crossing must be right at the Kauzmann point. We believe (speculate) 
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that the glass line should cross the "order - disorder" line right at the Kauzmann point if the experiments can be 
done accurately and the theory shall confirm it if the model can be solved exactly. It is of great interest to solve a 
solvable toy model to test this idea. The Kauzmann point is a point in which the magnetization and the entropy 
difference between solid and liquid phases changes sign. In this region the positive slope disorder dominated second 
peak segment joins the thermal fluctuations dominated negative slope melting segment. This is the region in which 
effects of the disorder and of the thermal fluctuation are roughly of the same strength. 

The replica symmetry breaking solution found here can be used to calculate the detailed properties inside the glass 
state. This however require generalization of the theory to include dynamics, since most of irreversible phenomena are 
time dependent. In particular it would be interesting to estimate the time scales associated with quenched disorder. 
This is left for a future work. Also we have considered only the three dimensional GL model here. It can be applied 
to superconductors with rather small anisotropy. It would be interesting to generalize the calculation to the Lawrence 
- Doniach model and to the two dimensional case describing thermal fluctuations and disorder in more anisotropic 
layered superconductors and thin films.. 
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Figure captions 

Fig. 1 

Schematic ax — r plane phase diagram of the vortex liquid with the |i/>| disorder only, ay is the LLL scaled 
temperature, while r is the \ip\ 2 disorder strength. The dotted line in the is the glass transition line. Below the 
line, the state is described by a replica diagonal matrix m a (, = mS a b, while above line the vortex state has a nonzero 
Anderson Edwards parameter. 



Fig. 



Phase diagram with the |?/>| 4 disorder only in the ax — q plane. The dotted line marks the replica symmetry 
breaking glass transition. The upper part above the " tricritical" point of this line is given by cq.(72), while the lower 
part below the "tricritical" point of the line is given by eq.(60) . The dashed line is also given by eq.(60), but does 
not correspond to a phase transition line. It is just a bifurcation line between two replica symmetric (the diagonal 
and the off diagonal) states. 

Fig. 3 

The H — T Phase diagram for the disorder only (q = 0) for three different disorder strength n = 0.12, 0.3 and 
0.08. The dashed dotted line is H C 2(T), the dotted line is the glass transition line, the solid line is the order-disorder 
(ODT or liquid -solid) transition line. Black blobs are the intersection points of the GT and the ODT lines. 

Fig. 4 

The schematic phase diagram of YBCO. Four distinct thermodynamic phases are: PL - pinned liquid (or vortex 
glass), PS - pinned solid (or the Bragg glass), S - solid (unpinned solid), and L - liquid (unpinned liquid). 

Fig. 5 

Comparison of the theoretical YBCO phase diagrams with the experiment. The dotted and the solid lines are 
the theoretical glass transition and the order-disorder transition lines respectively for the disorder strengths n = 0.12 
and q = 0. H m (T), H sp (T), H g (T) are the experimental the melting (ref.12), the second peak (ref. 3 3 ), and the glass 
transition (ref. 12) lines respectively. 
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